Two-dimensional imaging of the spin-orbit effective magnetic field 
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We report on spatially resolved measurements of the spin-orbit effective magnetic field in a 
GaAs/InGaAs quantum- well. Biased gate electrodes lead to an electric-field distribution in which 
the quantum-well electrons move according to the local orientation and magnitude of the electric 
field. This motion induces Rashba and Dresselhaus effective magnetic fields. The projection of the 
sum of these fields onto an external magnetic field is monitored locally by measuring the electron 
spin-precession frequency using time-resolved Faraday rotation. A comparison with simulations 
shows good agreement with the experimental data. 



In the reference frame of a moving electron, electric 
fields transform into magnetic fields, which interact with 
the electron spin and couple it to the electron's orbital 
motion, leading to spin-orbit (SO) interaction. In crys- 
tals lacking an inversion center such as GaAs, effective 
magnetic fields due to bulk inversion asymmetry (BIA) 
were predicted by Dresselhaus [l|. In heterostructures, 
structure inversion asymmetry (SIA) leads to an effec- 
tive magnetic field called Rashba term Both contri- 
butions have been studied extensively (for a review, see 
Ref. |3|) and are thought to play a crucial role in future 
spintronic devices, because the coupling of the orbital 
and the spin degrees of freedom opens a new way to spin 
manipulation, for example by flipping spins with oscil- 
lating electric fields [1, H H, 0] • The interplay between 
Dresselhaus and Rashba SO interaction has been stud- 
ied [1, H and proposed for use in a spin transistor [l(| • 
For finite electron wave numbers, SO interaction leads 
to a spin splitting at zero external magnetic field, which 
is observable as a beating of Shubnikov-de Haas oscil- 
lations [Ti|, Ei El EI El El E3 and can be used to 

experimentally determine the strength of the total SO 
interaction. The zero-field spin splitting also results in a 
spin-selective momentum scattering and can lead to spin- 
dependent photocurrents [l||, which allow the determi- 
nation of the ratio between the Rashba and Dresselhaus 
contributions by studying their directional dependence. 
Since SO interaction is the main reason for spin relax- 
ation in GaAs quantum well (QW) samples, the direc- 
tional dependence of spin relaxation reveals the relative 
strength of the Rashba and Dresselhaus terms [l||. A 
more direct method to get access to the SO fields is to 
impose a drift momentum on the conduction-band elec- 
trons by applying an in-plane electric field. This leads 
to two different effects: In a steady-state situation, the 
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spins are oriented along the SO field [2fJ, [21|, [23, l23| . If 
the spins of an ensemble of such drifting electrons are 
polarized by means of an optical pump pulse, a coher- 
ent spin precession about the SO field takes place, which 
has been observed using time-resolved Faraday rotation 
(TRFR) in a bulk GaAs epilayer 0]. Using TRFR, we 
have recently shown that the absolute values of both the 
Rashba and the Dresselhaus field can be determined in a 
two-dimensional electron gas by studying the precession 
frequency of electron spins as a function of the electron's 
direction of motion with respect to the crystal lattice [11 . 
An oscillating electric field was applied at an arbitrary 
direction in the plane of an InGaAs quantum well (QW) 
using two pairs of opposed electric gates arranged per- 
pendicularly to each other and enclosing a square area of 
QW electrons [see Fig. Ufa)]. 

Here, we study the spatial distribution of the total SO 
effective magnetic field between and outside the four gate 
electrodes that are used to generate the in-plane electric 
field. We show that the measured maps of spin-precession 
frequency can be well explained by assuming that the 
electrons move along the spatially varying electric field 
and that their spins perceive a SO effective magnetic field 
given by their local drift momentum. Depending on the 
orientation of an external magnetic field with respect to 
the crystal axis, the precession frequency becomes sen- 
sitive to drift momentum along a fixed in-plane direc- 
tion, allowing to spatially map the corresponding mo- 
mentum component. A simulation with fixed values for 
the Rashba and Dresselhaus coefficients as determined 
in Ref. [11 leads to good agreement with the measured 
precession-frequency maps. Specifically, the results con- 
firm that the electron drift momentum in the center of 
the gates can be described by the superposition of two 
perpendicular components given by the two orthogonal 
gate biases, as it was assumed in Ref. [!!■ Away from 
the center, however, the electric field and thus the spin- 
precession frequency start to vary. Large SO-induced 
effects on the spin-precession frequency can still be de- 
tected far away from the region between the gates, be- 
cause of the connection lines to the gates. This illustrates 
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the importance of a precise definition of the electrical 
connection scheme in a larger set-up and the difficulty of 
applying electric fields only at given positions of a two- 
dimensional electron gas. 

The samples are either a 20- or 43-nm-wide In- 
GaAs/GaAs QW as described in Ref. [25j (sample 1 
and 2). We experimentally determine B tot = |-B to t| = 
hu/g^B by measuring the electron spin- precession fre- 
quency v of optically excited conduction-band electrons 
using scanning TRFR (h denotes Planck's constant, g 
the electron g- factor, and ^b the Bohr magneton). With 
a first, circularly polarized pump pulse tuned to the 
absorption edge of the QW (A = 870 nm, average 
power 400 //W, pulse width of 2 ps, and repetition rate 
of 80 MHz), we create a spin polarization in the QW 
conduction band perpendicular to the QW plane. A 
pump-probe delay time At later, we probe the spin po- 
larization with a linearly polarized probe pulse (average 
power 60 //W) and monitor the rotation angle Of of its 
polarization plane. This Faraday rotation angle can be 
fit to 6 F = O exp (— Ar/TJ) cos (2twAt) to yield both 
the spin-precession frequency v and the spin-coherence 
time T$. By spatially moving the sample relative to the 
laser beams, the spin dynamics can be studied with a 
spatial resolution limited by the diameter of the beam 
in the focus, here approximately 15 /mi. Scanning Kerr 
microscopy has been used in Refs. [26j for spatially re- 
solving nuclear imprinting effects and in Refs [27t |28| for 
studying spin transport in GaAs epilayers. 

On top of the sample, four gates and the correspond- 
ing connection lines were defined using standard electron- 
beam lithography, evaporation of a 10-nm-thick Ti adhe- 
sion layer and a 80-nm-thick Au layer, and lift-off tech- 
niques. As shown in Fig. [IJa), opposite gate-pairs were 
connected to phase-locked oscillators, one each for the 
two perpendicular directions a and b. Two neighbor- 
ing gates were grounded, and by adjusting the voltage 
amplitudes V a and Vb on the other two gates, an oscil- 
lating electric field E(t) = E sin (2tt ft), f = 160 MHz, 
is induced in the QW plane. In the center of the four 
electrodes, Eq points at an angle, with respect to the a- 
axis, given by arctan Vb/V a . We have fabricated different 
samples in which the a-axis is oriented along either [110] 
or [100] of the semiconductor crystal. We define a co- 
ordinate system with x||[110] and j/||[110], and label the 
angle of E with the x-axis as tp, see Fig. QJb). As long 
as 1// is large compared to the mean collision time of 
the electrons in the QW (« 1 ps), the electrons adiabat- 
ically follow E(t) and their drift wave vector is given by 
k = m* /j,E(t)/hbar, with m* the effective electron mass 
and [i the electron mobility. Furthermore, the small colli- 
sion time prevents the electrons from leaving the laser fo- 
cus during the time 1//. Because of their non- vanishing 
A:- vector, the electrons are subject to SO effective mag- 
netic fields, which in two-dimensional systems are given 

by BEj 

B SIA = ^(M and B BIA = iqM, (!) 
g\x B \-k x J que \k x J 
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FIG. 1: (a) Microscope image of sample structure and 
schematic gate connections, (b) Relevant fields and angles. 
A magnetic field B cx t is applied at an angle 9 with respect to 
the [110] axis, and an electric field E at an angle tp. (c,d) Ori- 
entation of (c) Rashba and (d) Dresselhaus effective magnetic 
fields for the electron wave vector t on a unit circle. (e,f) 
Polar plot of A(6, tp) as a function of tp for 6 — 90° and 45°, 
as calculated for Bsia =-4mT and Bbia = 21mT (sample 1). 
Depending on 6, A is proportional to either the component of 
Eq perpendicular or approximately parallel to B C xt- 



with k = (k x , k y ) and a and (3 the Rashba and Dressel- 
haus coupling constants, respectively. The geometrical 
dependence on k is shown in Fig.[ljc) and (d). Whereas 
£?sia is always perpendicular to fc, Bbia points along 
±fe for the [100] and [010] directions. When analyzing 
the spin precession, both SO contributions can be added 
to an external magnetic field B ext [30l l3lj , which we ap- 
ply in the plane of the QW and at an angle 8 with the 
x-axis, as illustrated in Fig. [IJb). For all data presented, 
we choose B cxt = 0.958 T. 

If £>sia, -Bbia <§C -Bcxt, the total magnetic field B to t 
can be expressed by (2f| 

Btot(*) w B cxt + A{8, tp) sin {2k ft). (2) 

Here, t is the time delay between the electric-field oscilla- 
tion and the pump pulse (known up to an offset to, which 
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is constant in all experiments), and 

A{9, tp) = (B mA + B S ja) cos 9 sin tp 
+ (-Bbia - -Bsia) sin 9 cos tp. 

By probing v at different times t, we find B to t(t) to 
oscillate with t with the frequency /. Figure EJa) shows 
TRFR scans for two different t. This data is obtained on 
sample 1 in the center of the four gates that are aligned 
with the x- and y-axis. With V a = V& = 1 V, the electric 
field is oriented along [100] at the position of the laser 
focus. B cxt is applied along [110] (9 = 90°) . The TRFR 
signal monitors the coherent precession of the QW elec- 
tron spins about the total field B tot {t). Spins precess 
faster for t = -1.4 ns than for t= 1.7 ns. This variation of 
the precession frequency follows the oscillation of E{t), 
as becomes evident when plotting the fitted spin pre- 
cession frequency v as a function of t, see symbols in 
Fig. Eljb) . The data points fit to a sinusoidal oscillation 
at frequency / = 160 MHz (solid line). Except from a 
much weaker contribution at 2/ [25| . we do not observe 
higher harmonics in u{t), from which we conclude that 
the linear dependence of the effective spin-orbit magnetic 
field on k in Eq. [T] is valid. Specifically, terms cubic in k 
can be neglected. Such terms can be significant for the 
zero-field SO splitting [32| at the Fermi wave vector. In 
our case however, the effective SO magnetic field is deter- 
mined by the in-plane drift wave vector fc, which is much 
smaller than the quantized wave number perpendicular 
to the QW, and therefore negligible cubic contributions 
are expected. 

The obtained v{t) can be converted into B to t(t), the 
amplitude of its oscillation being given by A(9, tp) accord- 
ing to Eq. [21 By measuring A(9, tp) at varying and tp in 
the center of the four gates and comparing to Eq. (3J we 
determine Bsia = —4.2 mT, -Bbia = 21.6 mT for sample 
1, and Bsia = -8.5 mT, B B ia = 21.1 mT for sample 2. 
These values were obtained for a gate-modulation ampli- 
tude of 2 V, corresponding to E = 2900 V/m (2j|. 

In a geometrical interpretation, A{9, tp) is the projec- 
tion of Bsia + -Bbia onto B cxt . Because the SO fields 
depend linearly on Eq, A{9,tp) is proportional to the 
projection of the local field Eq onto an axis defined by 
(x,y) = ((Bbia - B S ta)sui0, (Bbia + B S ia)cos6»). At 
constant 9, a spatial map of A(9, tp) therefore directly im- 
ages one component of Eq. Figure [He) and (f) show the 
calculated A(9, tp) using the measured values for Bbia 
and Bsia as a function of tp for 9 = 90° and 45° for 
sample 1. In the first case, A is proportional to the 
projection of E onto the cc-axis, with an amplitude of 
Bbia — Bsia- For 9 = 45°, the projection is onto tp ~ 34°. 
For vanishing SIA, the highest visibility would be exactly 
&t<p = 45°. 

The electric field E is given by a superposition of the 
electric fields induced by the two pairs of electrodes along 
the a and b axis, and can be numerically determined with 
a partial differential-equation solver (e.g. pdetool in Mat- 
lab) with boundary conditions given by the voltages ap- 
plied to the gate electrodes and their connection lines. 
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FIG. 2: (a) Measured TRFR vs. At on sample 1 with an 
electric field oscillating at 160 MHz along the [100] axis, and 
6 = 90° . Two curves are shown for two different times t at 
which the electric field has opposite signs (open symbols). 
The solid lines are a fit to an exponentially decaying oscilla- 
tion. In (b), the obtained precession frequency v is plotted 
vs t (symbols). The solid line represents a fitted harmonic 
oscillation at 160 MHz, yielding the amplitude A(0,tp). 



Figure [3] shows a situation for sample 1 with the gate 
electrodes oriented along the x- and y-axis, and at two 
different bias configurations, where either the top and 
the right gate electrode [Fig. H^a)] or only the right gate 
electrode [Fig. [31(b)] were set to Vq = 1 V, and all other 
gates were grounded. In the center of the electrodes, Eq 
points along the [T00]-axis {tp — 225°) or the [Tl0]-axis 
{tp = 180°), respectively, as expected from simply super- 
posing two uniform fields along the x and y axis that are 
proportional to the two gate biases. 

At every point in the two-dimensional sample plane, 
the local electric field gives rise to a local fc-vector, which 
allows the calculation of A{9,tp) using Eqs. (P) and ([3]). 
The result of this simulation is shown in Fig. [3Jc) and 
(d). Considerable SO fields are also expected away from 
the four gate electrodes, because of the electric fields in- 
duced by the connection lines. The corresponding mea- 
surements are shown in Fig.[3je) and (f). The agreement 
with the simulation is good, except for the values of A 
measured close to a gate edge that are a factor of « 2 
lower than in the simulation. There, the simulation as- 
sumes perfect edges, leading to very high electric fields. 
In reality, the edges are rough and round, and the electric 
field is lower. Moreover, the simulations neglect the fact 
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FIG. 3: (Color online) (a,b) Simulations of the electric po- 
tential (color scale) and the electric field (arrows) for two dif- 
ferent gate potential configurations. (c,d) Simulated effective 
SO field and electric field (arrows). (e,f) Measured effective 
SO field on sample 1. 



that the gates are vertically offset from the QW by 30 nm, 
which compared to their lateral separation of 150 /im is, 
however, a negligible distance. 

As expected from Fig.QJe), electric fields along x lead 
to high SO contributions to v. In Figs. E^c) and (e), 
the connection lines along y induce strong electric fields 
along x in the upper left and the lower right quadrant of 
the sample. These fields lead to a fc-vector in x-direction, 
and, as visible from Fig. [He), to a high A. In contrast, 
close to the horizontal connection lines along x, the high 
electric fields in y-direction do not lead to visible SO 
fields, as in this situation, Bso is perpendicular to B cxt 
and therefore A = 0. 

For the situation with three grounded gates [Fig.|3ff)], 
outside of the four gates the electric field component 
along x is large in the full lower right corner and par- 
tially in the upper right corner, leading to a high value 
of A at these positions. Inside the gates, Eq is mainly 
aligned along — x and decays from right to left because of 
the shielding from the upper and lower gates which are 
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FIG. 4: (Color online) Measurement (left column) and simula- 
tion (right column) of the SO effective magnetic field between 
the four central gate electrodes (blue color) for sample 2. The 
left and the upper gate electrode are grounded, and biases V a 
and Vb are applied to the bottom and right electrode, respec- 
tively. The direction of the electric field in the center of the 
gate electrodes is indicated by black arrows. The calculated 
electric field is visualized as green arrows (right column). The 
magnetic field is applied along the [100] axis, which points to 
the right in the Figure. In the left column, the blue color is 
used for positions where no TRFR signal could be observed 
because the beams were blocked by the gate electrode. The 
color code is different than in Fig. [3] to visualize the sign of 
A. 
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at ground potential, which is directly observable in the 
measured and simulated maps of A. 

In Fig. IU we focus on the area between the four gate 
electrodes. We choose a configuration in which the gates 
and therefore the (a, 6)-axis are oriented at 45° to the 
(x, y)-axis. The measurements were taken on sample 2, 
and for the simulations the obtained values for £?bia and 
£?sia are used. We set = 45°, such that A is sensi- 
tive to the component of Eq approximately along B cxt , 
i.e. along the [100]-axis [see Fig. [TJf )] . The angle tp Q of 
the electric field in the center of the electrodes has been 
rotated by 30° in each step by varying the amplitudes 
V a = -V cos(45° + <fio) and V b = -y sin(45° + ip ) of 
the two phase-locked oscillators, connected to the bottom 
and the right gate electrode, respectively, with Vb = 1 V. 
This direction is indicated by an arrow between the mea- 
surement (left column) and the simulation (right col- 
umn). The simulation is again in good agreement with 
the measurement, except close to edges, where the sim- 
ulated electric field and therefore also the SO fields are 
higher than observed. 

In this configuration, A is sensitive to the horizontal 
component of Eq, and therefore, the highest A are mea- 
sured close to the right gate electrode (connected to Vb), 
provided that Vb is large enough. There, A is positive 
for <p = 225° and 195° (positive Vb) and negative for 
<p = 105° and 75° (negative Vb). As can be seen from 
the simulations in the right column of Fig [3l the sign of 
A correlates with the component of Eq along [100]. For 
ip = 135°, the left and right electrodes are grounded, and 
Eq points along the vertical direction in the center of the 
gates. On both the left and the right side of the center, 
Eq turns sideways towards the respective grounded lat- 



eral electrode, leading to positive and negative values for 
A on the two sides. Even though the bottom electrode 
induces strong electric fields, only a small A is observed 
close to the electrode, since the field is mainly oriented 
along [010], i.e. perpendicular to [100]. For <p = 165°, E 
has both positive and negative components along [100] 
close to the right gates, therefore leading to the appear- 
ance of both signs of A, as can be seen both in the simula- 
tion and in the measurement. In the corners of the square 
defined by the four gate electrodes, high A are measured 
as long as the two neighboring electrodes are on differ- 
ent potentials. There, the electric field is diagonal and 
therefore has substantial components along [100]. 

In conclusion, we have shown that the electric field 
induced by biased gate electrodes modifies the electron 
spin precession in an InGaAs QW through the presence 
of an SO effective magnetic field. Spatially resolved 
optical measurements of the spin precession frequency 
are in good agreement with numerically obtained spa- 
tial maps. They confirm that the change in spin pre- 
cession frequency is a measure of the projection of the 
drift momentum and thus the electric-field along an in- 
plane direction that is given by the Rashba and Dres- 
selhaus constants as well as the direction of B cxt - By 
adjusting the angle of B ext , different components of the 
electric field can be spatially mapped. Alternatively, if 
the electric field is known, this technique might allow to 
measure spatial variations of the Rashba and Dressclhaus 
constants [H, [34| . 
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